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Abstract. The possibility of intrinsic superconductivity in alkali-coated graphene monolayers has been 
recently suggested theoretically. Here, we derive the possible pairing symmetries of a carbon honeycomb 
lattice and discuss their phase diagram. We also evaluate the superconducting local density of states 
(LDOS) around an isolated impurity. This is directly related to scanning tunneling microscopy experiments, 
and may evidence the occurrence of unconventional superconductivity in graphene. 
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1 Introduction 

Graphene is a one-atom-thick layer of graphite [Ij . Due to 
its nearly ideal two-dimensional character and the linear, 
relativistic-like dispersion of its elementary excitations [2] , 
graphene is an intriguing material, where correlation and 
reduced dimensionality may conspire in favor or against 
various electronic instabilities. Indeed, it has been pro- 
posed that an electron liquid in a honeycomb lattice can 
be characterized by several ordered states, depending on 
doping and on the electron-electron interaction [3] . Among 
these competing orders, superconductivity (SC) could be 
stabilized by either topological disorder ^ or the prox- 
imity to an electronic topological transition (ETT) [5]. 
In this case, the symmetry of the underlying lattice may 
allow for an unconventional structure of the SC order pa- 
rameter, and possibly sizeable critical temperatures Tc, as 
is the case for the high- Tc cuprates [6l|7l[8]. 

The great interest in future graphene-based technol- 
ogy largely owes to the relatively high and easily tunable 
conductivity that can be realized in clean graphene sam- 
ples. Electronic correlations and the possible occurrence 
of ordered phases are therefore often neglected. Recently, 
a superconducting current has been observed to prop- 
agate through a superconductor-normal-superconductor 
(SNS) Josephson junction, where the N region consisted 
of a graphene layer [9,. This provided evidence of SC 
phase coherence in graphene single layers. Other carbon- 
based compounds are also known to sustain superconduc- 
tivity, sometimes with fairly large critical temperatures 
Tc- These include (a) highly oriented pyrolytic graphite 
(HOPG) [To], (b) the graphite intercalated compounds 
(GIC) [TTlfT^fTB] . which may be described as graphene 
sheets alternated by alkali layers, mainly acting as charge 



reservoirs, (c) quasi-one-dimensional carbon nanotubes, 
and (d) quasi-zero-dimensional alkali-doped fuUeres, or 
fuUerides ^14) . In particular, recent ab initio calculations 
support the idea that in-plane phonons may be respon- 
sible of SC in GIC [I l [T6 l [17 l [T l [l l [2g]. The quest for 
(quasi)intrinsic superconductivity in quasi-two-dimensional 
graphene is therefore well motivated, and several mecha- 
nisms have been proposed. These range from conventional, 
phonon-mediated superconductivity as in the GIC [TCI 
[TB1[T71[T^I191I201I21) . to unconventional, electronic mecha- 
nisms, including particularly the resonating valence bond 
(RVB) mechanism [SKH]- Moreover, it has been proposed 
[22] that SC may be realized in alkali-coated graphene sin- 
gle layers, where either electronic states belonging to the 
metallic bands get paired by means of graphene phonon 
modes, or graphene electrons couple via metallic plasmons 



Here, after classifying the symmetries of a SC order pa- 
rameter compatible with the honeycomb lattice of graphene, 
and discuss their stability and possible mixing as a func- 
tion of doping and coupling strengths, we study the local 
density of states (LDOS) of SC graphene around an iso- 
lated impurity [23] ■ In particular, we suggest that scan- 
ning tunneling microscopy (STM) measurements can dis- 
tinguish among the various available SC symmetries in 
graphene, thereby evidencing the possible unconventional 
nature of the SC order parameter, as is the case for the su- 
perconducting state of high-Tc cuprates 24, 2S], and pos- 
sibly of their unconventional normal state [2^12711^ . 
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2 Superconducting phase diagram 

We start by considering a model Hamiltonian H = Hq + 
Hi, where Hq = J2kx ^k\C\^\C]i\ describes the normal elec- 
tron liquid, with c^^ (cica) a creation (annihilation) oper- 
ator for a quasiparticle with wavevector k within the first 
Brillouin zone (IBZ) and band index A = ±, ^kA — E\^x~IJ- 
the tight-binding dispersion relation for band A, measured 
with respect to the chemical potential /i. Retaining hop- 
ping and overlap terms between nearest neighbor sites 
pg] . one has E]^x = At|7k|/(1 - As|7k|), where t = 2.8 eV 
and s = 0.07 are the nearest neighbor hopping and over- 
lap parameters, respectively [20], and 7k = Y^i=i e^^'^^ is 
the usual (complex) structure factor in momentum space. 
Here, 5i = a(l,V3)/2, 82 = a(l,-V3)/2, 5^ = a(-l,0) 
are the vectors connecting nearest neighbor sites in real 
space, where a = 0.142 nm is the C-C distance [2j. 

As for the pairing Hamiltonian Hi, we restrict to on- 
site (Von) and nearest-neighbor (y„„) interaction only. In 
terms of separate sets of creation and annihilation opera- 
tors a\^^, akcr (^kcr, ^kcr) for the A {B) sublattices, along 
the singlet channel and within the mean-field approxima- 
tion, one finds 

k 

+Ai{k){al^bl^^ - al^bl^^) + H.c, (1) 

where the two components of the order parameter are de- 
fined as 

k' k' 

-^i (k) = ^ ^ 7k-k' (&-k'iak't - b-k'tak'l) ■ (2b) 

k' 

While the onsite order parameter Aq is manifestly k- 
independent, the presence of an anisotropic intersite com- 
ponent ^i(k) opens the possibility of unconventional su- 
perconductivity. In particular, it is possible to factorize 
Z\i(k) in terms of the basis functions of the irreducible 
representations (irreps) of the point group Dq^ as 

2 

Z\i(k) =^zi„[0,(k)-f zV,(k)], (3) 

where 0o(k) -I- iV'o(k) — ^k/VS and 4>i{li) + i'0^(k) = 
^^k■5, _ ^ ;^^2,3) 311. Their subsets {0o}, {ipo}, 

{01,02}, {4'i,ip2} form a basis for the irreps Aig, Biu, 
E2g, and Eiu, respectively. One finds that the basis func- 
tions of both one-dimensional irreps vanish at the Dirac 
points K=§f(l,^), K' = |f(l,-^). Following stan- 
dard nomenclature]?], a nonvanishing component Aq, Aqi, 
and Aoe (£ = 1,2) will be termed s-wave, extended s- 
wave, and d-wave, respectively. We have analyzed the gap 
equations ([5]) as a function of temperature and chemical 
potential, as well as of the coupling constants. As is well 



known, when several competing symmetries are available 
for the order parameter, the gap equations factorize at 
Tc, and symmetry mixing is only possible below Tc |32j . 
Since the system is actually two-dimensional, the mean- 
field Tc only provides an upper bound for the Berezinskii- 
Kosterlitz-Thouless transition [33] . Fig. [1] shows the low- 
temperature mean-field phase diagrams for graphene, at 
various dopings, as a function of the coupling strengths. 
Solid lines separating filled regions are defined as the locus 
in the plane of reduced coupling constants, (Von/t, Vnn/t), 
such that the system exhibits SC with the relatively low 
critical temperature Tc = 4 K. The limiting case of pris- 
tine graphene (/i = 0, Fig. [Tl upper left panel) requires 
relatively high coupling strengths in order to develop SC 
even with a critical temperature as low as Tc = 4 K. 
At /i = 0, as T — >■ 0, one recover the phase diagram of 
Ref . [22] , where however only the s- and p + ip- waves sym- 
metries have been addressed. When /i ^ 0, one finds that 
Tc — >■ when Von, — >■ 0~, in agreement with results of 
Refs. [2^341135] . One then needs finite doping in order to 
have SC at moderately low coupling strengths (Fig. [TJ up- 
per right panel, /i = 0.5t). Upon doping, d-wave symmetry 
prevails at Tc close to the Van Hove singularities in both 
the valence and conduction bands, while s- or extended 
s-wave symmetry wins out close to the band edges. This 
is a generic effect of the proximity to an ETT, and has 
been described as a feature of high- Tc superconductivity 
in the cuprates 36 . At exactly the ETT within the con- 
duction band (Fig. [2 lower left panel, jjL = 1.07t), d-wave 
symmetry is maximally favored. It has been suggested [5] 
that fine-tuning the gate voltage in suspended graphene 
may achieve dopings of practical use. In particular, this 
may be even easier in uniaxially strained graphene, where 
the energy slope of the DOS increases with increasing 
strain [37]. Upon further increasing the chemical poten- 
tial (Fig. [Tj lower right panel, = 1.5t), the pure d-wave 
SC region shrinks, and eventually gets suppressed alto- 
gether. Assuming Von — (Z\o = 0), and neglecting band 
asymmetry (s = 0), one recovers the results of [ff'.T. In 
particular, Vnn — — 2.4i corresponds to Baskaran's pro- 
posal of resonating bond ordering in graphite [5]. 

3 Single impurity effects 

We next include the effect of a nonmagnetic localized im- 
purity, whose potential is given by 

where !^'ct(x) is a field operator at position x with spin 
projection a, and the energy [/q measures the impurity 
potential strength. In terms of the Nambu spinor ip\^ = 

('^kf ^kf "^-ki, ^-kj.^ , the imaginary time Green's func- 
tion is defined as tj(k, k',r) = — ((^k(T)</3k/ (0)), and obeys 
the Gor'kov matrix equation of motion 

^[(za;„ - /Ck)(5kq - V(k, q)]C/(q,k',iw„) =(5kk', (5) 
q 
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Fig. 1. (Color online) Low-temperature phase diagrams for 
the symmetry of the SC order parameter as a function of the 
coupling strengths. Left to right, top to bottom panels show 
the cases /i = (undoped), 0.5t, 1.07t (Van Hove singularity), 
and 1.5t. 



where fia;„ — {2n + 1)77^3^ is a fermionic Matsubara fre- 
quency. Here, /Ck is a 4 x 4 matrix associated with the 
mean-field SC Hamiltonian, 



K -^k 



whose 2x2 blocks are given by 



(£;k- - Sk+)/(27k) 
^i(-k) Ao 



i?k+) - M (i?k- - i?k+)/(27k) 



(6) 

(7a) 
(7b) 



Analogously, the impurity potential enters as a 4 x 4 block- 
diagonal matrix 

v(.,.',^C'"'o''''_„U,,)' 

whose 2x2 blocks have generic elements 

V^y{Kk') = a2t/o^L(x)^k'A'(x), (9) 

and ipk\{'x.) is the Bloch wavefunction for sublattice A = 
A, B employed in the tight-binding diagonalization of the 
pure sector of the Hamiltonian [35] ■ Eq. ^ is exactly sol- 
uble in the SC case, in the absence of impurities, while in 
the normal case and in the presence of localized impurities 
at highly-symmetric lattice positions it has been studied 
in [35] • In the general case, its solution can be expressed 
in terms of an appropriate T-matrix as 

5(k,k',ia;„) = (5kk'^/o(k, iw„) 

+5o(k, iw„)T(x; k, k', iw„)eo(k', lujn) (10) 

where t/o(k, ia;„) is the solution of the SC pure case. Fourier 
transforming back in real space via V'kA (r) and performing 



an analytic continuation to real frequencies, the local den- 
sity of states at position r can be expressed as p{r, lu) — 
— ifm G(r, r, cj). In what follows, we will use Gaussian 
pseudoatomic wavefunctions 0(r) = Zg cxp(— Z|r^/24a^)/(2aA/37r) 
{Zg = 11.2) [38] to expand the Bloch wavefunctions within 
the tight-binding approximation as ^kx{r) — N^^^^ (/)(r— 
j^A^jgjk Rj-^ where are vectors of the X = A,B sublat- 
tices. Placing the impurity on an A site, at r = say, and 
retaining only the zeroth order approximation, one has 
i'UAiO) ~ N-^/'^(j){0) and ■(/'kB(O) « 0. In this limit, the 
impurity potential matrix simplifies to V(k, k') = Vo = 
(l)'^{0)a'^UQN~^a3 (g) (ctq + o'3)/2, where ao, as are Pauh 
matrices. Correspondingly, the analytically continued T- 
matrix assumes the k-independent form T(r = 0;uj) = 
VoA^^^(0, cj), where zeroes of the determinant of 7M(0, w) = 
1 — X]q^o(q, w)Vo as a function of w are connected with 
impurity-induced bound states or well-defined resonances. 

Fig. [2] shows the LDOS on an atomic site (r = 0) for 
pure SC graphene {Uq = 0). Here, and in what follows, 
we shall consider slightly doped graphene (/i — 0.5t), in 
order to stabilize the SC phase without requiring exceed- 
ingly high coupling strengths. This corresponds to hav- 
ing the Fermi surface in the shape of two disconnected 
rings, centered around either inequivalent Dirac point. In 
the s-wave case (Fig. [51 top panel), the LDOS assumes a 
typical BCS-like shape, with well-defined coherence peaks 
located at w « ±Z\o/2, and no spectral weight is available 
within the gap. Away from the gapped region, the LDOS 
is not affected by superconductivity, and in particular the 
Van Hove singularities remain intact. Qualitatively simi- 
lar results have been obtained in the case of an extended 
s-wave order parameter. This is because, despite the non- 
trivial k-dependence, the order parameter vanishes pre- 
cisely at the Dirac points, viz. where the bands vanish if 
fi — 0. On the other hand, in the d-wave case (Fig. |21 bot- 
tom panel), the onsite LDOS vanishes linearly at w = 0, 
and the coherence peaks are considerably reduced. This is 
typical of unconventional superconductors, and is in fact 
a hallmark of d-wave superconductivity in the cuprates 
[27]. In Fig. [2] we have assumed in-phase d-wave com- 
ponents for the order parameter, i.e. An = A12. Other 
phase relations are possible, such as A12 — Aiie~^^^^. In 
the latter case, ^i(k) vanishes only at the Dirac points. 
Therefore, in the case of doped graphene, the SC order 
parameter has no node along the Fermi surface, and one 
finds a fully gapped LDOS, qualitatively similar to the 
s-wave case. This is in agreement with the results of [7]. 
However, while an out-of-phase order parameter is numer- 
ically favored in the case of pure d-wave symmetry, we 
find that a nonzero s-wave component stabilizes in-phase 
d-wave superconductivity. 

We consider now the effect of a single, localized, non- 
magnetic impurity, located on a lattice site. In the normal 
state, we recover the results of [35] for a site-like impurity. 
Fig. |3] shows our results for the LDOS in the SC state. In 
the s-wave case (Fig. [3] top panel), the onsite LDOS is 
characterized by the opening of a gap and the formation 
of well-defined coherence peaks, while the Van Hove sin- 
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Fig. 2. (Color online) LDOS for pure SC (solid line) and 
normal-state graphene (dotted line). Top panel refers to s-wave 
symmetry, while bottom panel refers to d-wave symmetry, with 
in-phase components of the order parameter. In both cases, 
p = 0.5t and Uo = 0. 

gularity gets suppressed, as in the normal state. On the 
other hand, the LDOS on a nearest-neighbor site gets en- 
hanced with respect to its normal state counterpart, and 
the Van Hove singularities undergo no significant reduc- 
tion. More importantly, by analyzing the determinant of 
the matrix A4{r, uj), one finds that no bound state is pos- 
sible within the gapped region in the s-wave case. In the 
d-wave case (Fig. [31 bottom panel), while a similar analysis 
applies as in the s-wave case at large energies, the situa- 
tion is completely different within the gapped region. One 
finds that sufficiently high impurity coupling strengths 
{—3.8t ^ Uq < —0.8t, in the case under study) allow 
the formation of well-defined bound state pairs within the 
gapped region. (That bound states should come in pairs 
follows from the particle-hole mixing, characteristic of the 
SC state [17].) This is indeed apparent from Fig. [3] (bottom 
panel), with the appearance of sharp peaks in the onsite 
LDOS. Such a result is similar to the d-wave SC state of 
the cuprates, and has indeed been verified experimentally 
in STM experiments in Bi2212 [24l[25] . 

4 Conclusions 

In summary, we derived the possible symmetries of a SC 
order parameter compatible with the honeycomb lattice 
structure of graphene. We discussed the SC phase diagram 
at low-temperature as a function of doping and interaction 
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Fig. 3. (Color online) LDOS in the presence of a localized 
impurity for SC (solid lines) and normal-state graphene (dotted 
lines). Green lines refer to LDOS on top of the impurity (x = 
0), while blue lines refer to LDOS on a nearest-neighbor site 
(x = Si). Red lines refer to the on-site LDOS (x = 0) in 
the SC phase, in the absence of impurities. Top panel refers 
to extended s-wave symmetry, while bottom panel refers to 
d-wave symmetry. In both cases, jj, = 0.5t and Uo ~ — 0.75t. 

strengths. In particular, we showed that unconventional 
pairing may stabilize close to an ETT, with possible ad- 
mixtures of a subdominant s-wave contributions. 

We then evaluated the LDOS around an isolated im- 
purity in the SC phase, and suggested that STM experi- 
ments may detect the occurrence of unconventional pair- 
ing in intrinsically or proximity-induced superconducting 
graphene. This will be helpful in determining the pairing 
mechanism of superconducting graphene. We have lim- 
ited our study to isolated nonmagnetic impurities. Such 
a condition may be realized in the limit of low impurity 
concentration, where it is safe to neglect impurity effects 
on the overall transport properties of graphene. 
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